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(Dated: March 28, 2012) 

We measured separate cosmic-ray electron and positron spectra with the Fermi Large Area Tele- 
scope. Because the instrument does not have an onboard magnet, we distinguish the two species 
by exploiting the Earth's shadow, which is offset in opposite directions for opposite charges due 
to the Earth's magnetic field. We estimate and subtract the cosmic-ray proton background using 
two different methods that produce consistent results. We report the electron-only spectrum, the 
positron-only spectrum, and the positron fraction between 20 GeV and 200 GeV. We confirm that 
the fraction rises with energy in the 20-100 GeV range. The three new spectral points between 100 
and 200 GeV are consistent with a fraction that is continuing to rise with energy. 



PACS numbers: 96.50.S-, 95.35. +d, 95.85.Ry, 98.70.Sa 

Introduction. Positrons were discovered in cosmic rays 
(CRs) in 1964 [1]. Subsequent measurements through the 
1960s and 1970s [2, 3] showed that the positron fraction, 
e+/(e+ + e~), decreases with energy up to ~10 GeV. 
See [4] and references therein for a summary of these 



early experiments. Balloon payloads through the 1980s 
measured CR electrons and positrons with inconsistent 
results [5, and references therein]. Starting in the mid- 
1990s, HEAT [6-10], CAPRICE [11, 12], and TS93 [13] 
measured the CR positron spectrum and fraction up to ~ 



50 GeV, indicating a predominantly decreasing positron 
fraction with increasing energy. However, a small excess 
in the positron fraction above ~ 7 GeV was detected by 
HEAT and also seen in CAPRICE data, as well as by 
AMS-01 [14]. Recently, the PAMELA instrument has 
measured a positron fraction that increases with energy 
above ~10 GeV [15, 16] with high precision, confirming 
the indications seen in the earlier data. 

The best established mechanism for producing CR 
positrons is secondary production: CR nuclei interact 
inclastically with interstellar gas, producing charged pi- 
ons that decay to positrons, electrons, and neutrinos. 
However, this process results in a positron fraction that 
decreases with energy [4, 17]. The origin of the rising 
positron fraction at high energy is unknown and has been 
ascribed to a variety of mechanisms including pulsars, 
CRs interacting with giant molecular clouds, and dark 
matter. See [18, 19] for recent reviews. 

The Large Area Telescope (LAT) is a pair-conversion 
gamma-ray telescope onboard the Fermi Gamma-ray 
Space Telescope satellite. It has been used to measure 
the combined CR electron and positron spectrum from 
7 GeV to 1 TeV [20, 21]. The LAT docs not have a mag- 
net for charge separation. However, as pioneered by [22] 
and [23], the geomagnetic field can also be used to sepa- 
rate the two species without an onboard magnet. Miillcr 
and Tang [23] used the difference in geomagnetic cut- 
off for positrons and electrons from the east and west 
to determine the positron fraction between 10 GeV and 
20 GeV. As reported below, we used the shadow im- 
posed by the Earth and its offset direction for electrons 
and positrons due to the geomagnetic field, to separately 
measure the spectra of CR electrons and positrons from 
20 GeV to 200 GeV. In this energy range, the 68% con- 
tainment radius of the LAT point-spread function is 0.1° 
or better and the energy resolution is 8% or better. 

Region selection and exposure calculation. The Earth's 
magnetic field significantly affects the CR distribution in 
near-Earth space. At energies below ~10 GeV, a signifi- 
cant fraction of the incoming particles are deflected back 
to interplanetary space by the magnetic field ("geomag- 
netic cutoff"). The exact value of the geomagnetic cutoff 
rigidity depends on the detector position and viewing 
angle. In addition to the geomagnetic cutoff effect, the 
Earth blocks trajectories for particles of certain rigidities 
and directions while allowing other trajectories. This re- 
sults in a different rate of CRs from the east than the 
west (the "east-west effect") [24-26]. 

Figure 1 shows example trajectories for electrons and 
positrons. Positive charges propagating toward the east 
are curved outward, while negative charges are curved 
inward toward the Earth (Figure 1). This results in a 
region of particle directions from which positrons can ar- 
rive, while electrons are blocked by the Earth. At each 
particle rigidity there is a region to the west from which 
positrons are allowed and electrons are forbidden. There 
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FIG. 1: Examples of calculated electron (red) and positron 
(blue) trajectories arriving at the detector, for 28 GeV parti- 
cles arriving within the Equatorial plane (viewed from the 
North pole). Forbidden trajectories are solid and allowed 
trajectories are dashed. Inset: the three selection regions 
(electron-only, positron-only, and both-allowed) for the same 
particle energy and spacecraft position as the trajectory traces 
(viewed from the instrument position in the Equatorial plane) . 



is a corresponding region to the east from which electrons 
are allowed and positrons are forbidden. The precise size 
and shape of these regions depend on the particle rigidity 
and instrument location. 

We used a high-precision geomagnetic field model (the 
2010 epoch of the 11th version of the International Ge- 
omagnetic Reference Field [27]) and a publicly available 
code [28] to trace charged particle trajectories in the mag- 
netic field and determine allowed vs. forbidden regions 
for each species. We previously used the same magnetic 
field model and tracer code to perform a precise compar- 
ison between predicted and measured geomagnetic cutoff 
rigidities for the Fermi LAT orbit, finding that the tracer 
code accurately predicts the geographical distribution of 
the geomagnetic cutoff [29]. We used a static 2010 model 
for all of the data analyzed here, which were recorded be- 
tween June 2008 through April 2011. 

Each particle trajectory is traced backward from the 
spacecraft until it reaches 20 Earth radii from the Earth 
center or reaches the Earth's atmosphere, which we ap- 
proximate with a 60 km thickness (Figure 1). If the tra- 
jectory reaches 20 Earth radii, it is an allowed trajec- 
tory. If it reaches the atmosphere, it is a forbidden tra- 
jectory. We calculate electron-only, positron-only, and 
both-allowed (control) regions for each 30 s time step us- 
ing the instantaneous spacecraft latitude and longitude 
and the nominal orbital altitude of 565 km. The regions 
are determined for each energy bin, with 10 logarithmi- 
cally spaced energy bins spanning 20-200 GeV. The 30 s 
time step (in which the spacecraft travels ~ 2° longitude) 
is sufficient to achieve a finely sampled distribution of 



instantaneous regions and exposures. Although we use 
binned position data for the exposure calculation, we use 
the instantaneous spacecraft position at the time of each 
event to determine which region it lies in, so the event 
selection is not affected by the finite step size. 

We define the "deflected horizon" (inset of Figure 1) 
to be the curve that separates the allowed and forbidden 
regions for a particular energy, charge, and spacecraft 
position. The curve represents the position of the Earth 
horizon as "deflected" by the geomagnetic field, with de- 
flection occurring in one direction for electrons and in the 
opposite direction for positrons. At 20 GeV the curves 
differ significantly from the actual Earth horizon. At 
higher energies, the curves asymptotically approach the 
undcflcctcd horizon. The positron-only region is defined 
to be the region above the positron horizon and below the 
electron horizon. A corresponding definition is used for 
the electron-only region. The region above both instan- 
taneous horizon curves and below a nadir angle of 130° is 
the control region, where both species are allowed. The 
control region is truncated at 130° because the additional 
statistics gained using a larger region are unnecessary for 
this analysis and require more processing time. The re- 
gions vary with spacecraft position and particle energy. 

We determine the integrated exposure (cm 2 sr s) for 
each of the three regions by integrating the energy- and 
direction-dependent effective area of the instrument over 
the solid angle of each region and over the livetime of the 
data set. The numerical integration is performed with 
0.5° pixels. For each particle, we use the reconstructed 
energy and actual spacecraft position at the time of the 
event to determine the deflected electron and positron 
horizon curves. We then calculate the nadir angle dif- 
ference between the particle arrival direction and the de- 
flected horizon: D± = N — N±, where N is the nadir 
angle from which the event arrived and N± is the nadir 
angle of the positron (electron) horizon, evaluated at the 
azimuth direction from which the event arrived. As de- 
scribed below, we use D± to select one set of CR elec- 
trons and another set of CR positrons and to remove 
atmospheric electrons and positrons. 

When a CR nucleus interacts in the Earth's atmo- 
sphere, it produces a shower of secondary particles in- 
cluding gamma rays, electrons, and positrons. At the 
high energies considered here, the secondaries are closely 
collimated along the direction of the primary. The limb 
of the atmosphere is therefore visible in gamma rays as a 
bright, thin (1-2° wide) ring of emission. [30]. The peak 
emission is produced at a nadir angle of 68°, which is 
1.3° (^60 km) above the Earth surface. 

Atmospheric electrons and positrons are produced 
with angular distributions similar to that of atmospheric 
gamma rays. Unlike the gamma rays, they are deflected 
by the geomagnetic field before reaching the LAT and 
arrive from the direction of the deflected horizon curves 
rather than the actual Earth horizon. After accounting 
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FIG. 2: Distribution of D+, for events with D- < -2°. The 
spike at D+ = 0° is due to atmospheric positrons. The events 
with D+ > 0° are CR positrons plus residual proton back- 
ground events. The events with D+ < 0° are residual back- 
ground events. 

for the deflection, however, they are also detected as a 1- 
2° wide ring of emission (Figure 2). Like the gamma-ray 
emission, the peak atmospheric positron emission seen by 
the LAT is produced at an impact parameter of ^60 km 
above the Earth surface. The detection of this popula- 
tion, and its narrow distribution at D + = 0, confirms 
that the magnetic field model and particle tracer code 
are accurate. As a cross check, we also ran the trajectory 
tracer code with the 1995 magnetic field model. Chang- 
ing from the 2010 to the 1995 model, D± changes by less 
than 1° for most events, and the distribution shown in 
Figure 2 changes negligibly. 

To remove atmospheric electrons and positrons, which 
are concentrated near D + = and £>_ = 0, we trim the 
edges of our selection regions by a "trim width" W in 
nadir angle. We evaluated the residual contamination by 
atmospheric electrons and positrons by performing the 
full positron flux measurement as a function of W . The 
contamination decreases significantly from W = 0° to 
W = 2° and is smaller than statistical fluctuations for 
W > 4°. The contamination is <10% for W = 2° and 
<3% for W = 4°. Below 100 GeV we trim by 4°. Above 
100 GeV we trim by 2°; the systematic effect from resid- 
ual secondaries with this less conservative trimming is 
larger (10%) but still smaller than the statistical uncer- 
tainty in this energy range. As explained below, we have 
included these estimates of the residual atmospheric con- 
tamination in the systematic uncertainties. This effect is 
included for the positron flux but is negligible for the 
larger electron flux: negligible variation with W is seen 
in the electron flux. 

Run and event selection. Data collected in the LAT 
standard mode of operation (survey mode) are not suit- 



able for this analysis. In survey mode, the LAT points 
away from the Earth, keeping the electron-only and 
positron-only regions (which are centered on the Earth 
horizon) outside of or at the edge of the field of view. 
However, the LAT has observed in special modes with the 
horizon in the field of view. Direct limb-pointed obser- 
vations have been performed for instrument commission- 
ing and calibration [30]. Pointed observations have also 
been performed in which a particular celestial direction 
is tracked for several hours or several days. During these 
observations the target is tracked until it is occulted by 
the Earth, during which time the spacecraft tracks along 
the limb of the Earth, with a 50° offset toward zenith, 
until the target re-emerges from occultation. We used 
775 runs from these observation categories with a total 
observation time of 39.0 days. A LAT run is an interval 
of data taking that typically spans one orbit. 

The recorded events are dominated by a background 
of CR protons and heavier nuclei. Here we used the 
same particle classification technique that we developed 
to measure the combined e + + spectrum [20, 21] 
to separate CR leptons from this background. We use 
events in the 20-200 GcV range. Below 20 GcV the 
electron/positron sensitivity is small due to the event 
filter applied onboard the LAT: any event that deposits 
20 GeV in the calorimeter passes this filter, but below this 
energy charged particles are significantly suppressed in 
favor of gamma rays [21]. Above 200 GeV, the electron- 
only and positron-only regions become very small; this 
combined with the steeply falling spectra makes the num- 
ber of detected events prohibitively small. 

Background subtraction. The CR hadron flux (pre- 
dominantly protons) between 20 GeV and 200 GcV is 
2-3 (depending on energy) orders of magnitude larger 
than the electron flux and 3-4 orders of magnitude larger 
than the positron flux [31]. Our lepton/hadron separa- 
tion criteria suppress the hadron background to a level 
that is between 4% and 20% (depending on energy) of 
the combined electron+positron flux [20, 21]. 

We used two independent methods to estimate the 
residual proton contamination remaining after the selec- 
tion cuts described above. The first background estima- 
tion method uses fits applied to the flight data alone; it 
does not require simulating the detector. The transverse 
size of showers in the calorimeter provides significant dis- 
crimination power for separating electromagnetic show- 
ers from hadronic interactions. In this method we use 
slightly relaxed selection criteria to increase the statistics 
of both leptons and hadrons and fit the distribution of 
transverse shower size to statistically estimate the num- 
ber of signal (and background) events present in the three 
regions. 

Distributions of the transverse shower size are pro- 
duced for each energy bin and for each region. They 
are fit with a function that represents the sum of signal 
and background (Figure 3). The number of leptons in the 
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FIG. 3: Transverse shower size distribution in the electron- 
only region. In the positron-only region, the number of events 
with small transverse shower size is smaller, but the mean and 
width of the distribution are similar. See [32] for additional 
figures. 

sample is determined by integrating the signal function 
up to 35 mm, where we truncate the integral to reduce 
the effect of the distribution tail. The signal and back- 
ground shower size distributions are each parameterized 
with a Gaussian. To evaluate the uncertainty due to 
small differences between real data distributions and the 
Gaussian approximation, we applied this fit-based tech- 
nique to a sample of events collected in survey mode to 
determine the e + + e~ spectrum, and compared it to the 
previously published LAT measurements [20, 21]. The 
spectrum determined with the fit-based background sub- 
traction method agrees with the published one within 
5%. Tests on Monte Carlo simulations confirm that the 
choice of Gaussian functions is sufficient to obtain the 
number of signal events at the few-percent level. 

The fit is well constrained in the control (both species 
allowed) region where the number of events is large and 
the signal-to-background ratio allows us to easily distin- 
guish protons from electrons/positrons. In the other two 
regions (electron-only and positron-only) the statistics 
decrease rapidly with energy and the signal peak pa- 
rameters are not well constrained. This is particularly 
important in the positron region where the background 
distribution has a larger overlap with the signal. How- 
ever, the distribution of shower size depends primarily 
on particle type (lepton or hadron) and on the angle 9 
between the particle track and instrument axis. The elec- 
tron and positron regions correspond to particular 9 dis- 
tributions (the two distributions are very similar to one 
another) which we reproduce by assembling a reference 
set of events from the control region with the same 9 dis- 
tribution. We then fix the mean and sigma parameters 
of the signal Gaussian to the values obtained by fitting 



this reference data set and fit only the normalization to 
the data sets from the positron-only and electron-only 
regions. No assumptions are made on incoming angle de- 
pendence of the background Gaussian, whose parameters 
are always free. 

We evaluated the uncertainty introduced by this fit 
stabilization procedure by changing the mean and sigma 
parameters (within ±3 times the uncertainty of the fit 
values from the reference data sets) and calculating the 
resulting variation in the number of signal events. The 
maximum deviation of the number of signal events is 4% 
and is included in the systematic uncertainty for each 
energy bin. 

The second method of estimating the residual proton 
background is based on a Monte Carlo (MC) simulation. 
The LAT detector MC is based on GEANT 4 [33] and 
has been validated using comparisons to flight data and 
pre-flight beam test data. We produced a set of 40 bil- 
lion Monte Carlo protons between 4 GeV and 10 TeV 
(~3% of which pass the trigger and onboard filter) with 
a power- law distribution in energy with an index of 1.5. 
We also simulated 150 million electron events between 
4 GeV and 1 TeV (-13% of which pass the onboard fil- 
ter) with a spectral index of 1. The spacecraft orbit in the 
simulation matched the actual flight orbit for the data set 
used. We re-weight the MC proton flux to the spectrum 
measured by the AMS-01 experiment (index 2.78) [34], 
and the MC electron flux to the spectrum measured by 
the Fermi LAT [20, 21] (index 3.08). We apply the same 
analysis cuts to the MC and flight data in order to esti- 
mate the residual proton contamination. Roughly 0.1% 
of the protons that pass the onboard filter also pass the 
analysis cuts we designed to reject them. We traced the 
trajectories of surviving events and removed events with 
trajectories blocked by the Earth. The accuracy of the 
MC proton rate has been validated by comparing a vari- 
ety of distributions between MC and flight data. In par- 
ticular, we inverted individual cuts in order to produce 
samples with an enriched proton background contribu- 
tion. The rate of MC protons agrees with flight data 
within ~ 8%. 

Results. The common systematic uncertainty in the 
CR electron/positron flux for the two background sub- 
traction methods arc: ±5% effective area; ±5% on- 
board filter efficiency in the two lowest energy bins (20.0- 
31.7 GeV); and +°| below 100 GeV, +^ % above 100 GeV 
due to atmospheric lepton contamination. The system- 
atic uncertainty of atmospheric positron contamination 
is asymmetric because such contamination can only cause 
us to overestimate, not underestimate, the CR positron 
flux. Systematic uncertainties present only for the flight 
data fitting method are 5-10% (depending on energy and 
region) due to discrepancies between the fit shape and 
the actual distribution and 2-4% due to using the ref- 
erence 9 distribution. Components present only for the 
MC method arc —8% due to discrepancies between MC 



protons and flight data, and 2-10% CR proton spectral 
index uncertainty, evaluated using three recent CR pro- 
ton spectrum measurements (BESS [35], AMS-01 [34], 
and PAMELA [36]). Summing systematic uncertainty 
components in quadrature (we expect them to be un- 
correlated), we estimate 10-16% (depending on energy) 
uncertainty in the positron flux for the fit method and 
8-19% for the MC method. To determine the final un- 
certainty of each spectral point we add the statistical un- 
certainty in quadrature; for the MC method this includes 
both signal and background statistics. 

Our best estimates of the spectra are shown in Fig- 
ure 4 and Table I. We chose the fit method for all energy 
bins except the highest, because this method has slightly 
smaller uncertainty. For the highest energy bin we use 
the MC method because the statistics are not sufficient 
for fitting. Results from the MC method are consistent 
within the total (systematic + statistical) uncertainty of 
the fit method in all of the bins in the control and electron 
regions and in seven out of nine bins in the positron re- 
gion. The maximum single-bin difference is 1.5 times the 
total uncertainty. The positron fraction ( j( c +) + j| c -) ) is 
shown in Figure 5. Error propagation is performed inde- 
pendently for the fraction and for the individual spectra; 
the contribution due to effective area uncertainty cancels 
in the fraction. The small difference between the com- 
bined e + + e~ flux we measured in the control region 
and the total flux we reported previously is due to in- 
strument response functions that have been updated to 
account for "ghost events" [37]. 



Energy J(e+) x 10 5 J(e~) x 10 4 
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TABLE I: Flux (GeV -1 rrT 2 s _1 sr _1 ) and positron fraction 
as a function of energy (GeV). Uncertainties are istat isys. 

The spectrum measured in each of the three regions be- 
tween 20 and 200 GeV is well described by a power law. 
The fit to the positron spectrum is (2.02i0.22 x 10~ 3 



GeV" 



sr 



)-2.77±0.14 whUe the fit tQ 



1. 20GoV > 

the electron spectrum is (2.07i0.13 x 10~ 2 GeV -1 m~ 2 



sr 



20GoV 



) 319 ± 07 . The uncertainties are deter- 



mined by including the total (statistical plus systematic) 
uncertainty of each energy bin. The fitted indices are con- 
sistent with the index we reported previously for the total 
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We find that the positron fraction increases with en- 
ergy between 20 and 200 GeV, consistent with results 
reported by PAMELA [15, 16]. Future measurements 
with greater sensitivity and energy reach, such as those 
by AMS-02 [40] , are necessary to distinguish between the 
many possible explanations of this increase. 
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in Japan, and the K. A. Wallenberg Foundation, the 
Swedish Research Council and the National Space Board 
in Sweden. Additional support from INAF in Italy and 
CNES in France for science analysis during the opera- 
tions phase is also gratefully acknowledged. 



FIG. 4: Energy spectra for e + , e~, and e + + (control re- 
gion). In the control region where both species are allowed, 
this analysis reproduces the Fermi LAT results reported pre- 
viously for the total electron plus positron spectrum [20, 21] 
(gray). Previous results form HEAT [9] and PAMELA [38] 
are shown for reference. The bottom panel shows that the 
ratio between the sum and the control flux is consistent with 
1 as expected. 
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FIG. 5: Positron fraction measured by the Fermi LAT and 
by other experiments [7, 14, 16]. The Fermi statistical uncer- 
tainty is shown with error bars and the total (statistical plus 
systematic uncertainty) is shown as a shaded band. 



electron plus positron spectrum (3.08±0.05) [20, 21]. 

Conclusion. We measured the CR positron and elec- 
tron spectra separately between 20 and 200 GeV, using 
a novel separation technique which exploits the charge- 
dependent displacement of the Earth's shadow due to 
the geomagnetic field. While the positron fraction has 
been measured previously up to 100 GeV [15, 16] and 
the absolute flux has been measured previously up to 
50 GeV [9, 39], this is the first time that the absolute CR 
positron spectrum has been measured above 50 GeV and 
that the fraction has been determined above 100 GeV. 
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MEASUREMENT OF SEPARATE COSMIC-RAY 
ELECTRON AND POSITRON SPECTRA WITH 

THE FERMI LARGE AREA TELESCOPE. 

SUPPLEMENTARY ONLINE MATERIAL. 

Shower transverse size distribution 

The distribution of the transverse size of showers in 
the calorimeter is used to determine the number of sig- 
nal events in the three regions used in this analysis: the 
electron-only, the positron-only, and the control region 
(where both species arc permitted) . In the paper we used 
the electron-only case between 31.7 and 39.9 GeV to ex- 
plain the technique. Here we show also the positron-only 
case in Figure 6. Figure 7 shows another example in the 
energy interval from 126.2 GeV to 158.9 GeV, which is 
the highest energy bin in which the fitting technique is 
applied. The distributions of events in the control region 
for both energy intervals are shown in Figure 8. 




FIG. 6: Transverse shower size distribution in the electron-only region (left panel) and positron-only region (right panel) for 
the energy interval from 31.7 GeV to 39.9 GeV. 
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FIG. 7: Transverse shower size distribution in the electron-only region (left panel) and positron-only region (right panel) for 
the energy interval from 126.2 GeV to 158.9 GeV. 




FIG. 8: Transverse shower size distribution in the control region in the two energy intervals 31.7 - 39.9 GeV (left panel) and 
126.2 - 158.9 GeV (right panel). 



